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A  Simple  Real-time  Model 
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Important  features: 

•  Nonlinear  System 

-  Hull  orientation 

-  Road-arm  rotations 

Ground  Contact 

-  Becker  law 

-  Constant  SA 

Jounce  Bumper 

-  Exponential  stiffness 


Tracked  Vehicle 
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Simplifying  assumptions: 

•  Track  loop  “ignored” 

-  Road  wheels  rotate  in  unison 

-  Specified  L/R  track  speeds 

•  Ground  contact 

-  “Spherical”  wheels 

-  Single  point  look-up 


Tracked  Vehicle 


•  smooth  terrain 

•  curvature1  »  RW  radius. 
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Force  Elements: 
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Terrain: 

•  Scalable  bump 

•  Height,  wavelength 

Maneuver: 

Track  speeds 
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Case  1 : 

•  Model  integrated  with 
error  control  (AS)  and  at 
200  Hz,  500Hz,  and 
1kHz 

•  Bump  height  and 
wavelength  <1 , 5>m 

•  Trajectories  are 
indistinguishable. 
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Results  :  Case  1 


Road  Arm  Rates  4  Road  Arm  Moments 

x  10 
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•  Peak  Loads 

—  AS  -1.31x105Nm  xio4  Road  Arm  Moments 


-  200Hz  -1 .21  xIO5  Nm 

-  500Hz  -1 .09x1 05  Nm 

-  1kHz  -1 .25x1 05  Nm  « 
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•  Results  vary  chaotically  | 

-  8%  1 

-  17% 

-  5% 


•  Aliasing  also  contributes. 
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Case  2: 

•  Again  model  integrated 
with  error  control  (AS) 
and  at  200  Hz,  500Hz, 
and  1  kHz 

•  Bump  height  and 
wavelength  <1 , 2>m 

•  Initial  terrain 
assumption  violated! 

•  Results  are  unstable. 
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Numerical  Temporal  Integration 


Runge-Kutta  family  of  explicit  integrators 

-  Forward  Euler,  RK3,  RK4,  ode23,  ode45,  Dormand-Prince,  etc. 

Basics 

-  Approximate:  y(t0+h)  where  y’  =  f(t,y) 

-  Taylor  Series  Expansion: 

y(t0  +  h)  —  y(t0 )  +  hy'(t0 )  +  -y  y"(t0 )  +  0(h3 ) 
y(t0  +  h)  =  y(t0)  +  hf  +  *Y^  +  0(h3) 

4f_  _  dt_  _  f  f  ,  f 

dt  dy  dt  dt  dt  J  yJ  J  t 

y(t0+h)  —  y(t0)  +  hf  +  -y(/r  +  fyf)  +  0(h3) 
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1  Numerical  Temporal  Integration 

MODELING  rho  SIMULATION,  TESTING™  VALIDATION  ’ 


-  Taylor  Series  Expansion: 

y(t0  +  h)  =  y(t0)  +  hf  +Iy(ft  +  fyf)  +  0(h3) 

-  Function  Evaluations: 

K=f(to^yo)  =  f 

k2=  f  (?o  +  ha2,y0  +  hb2lkx)  =  f  +  ha2ft  +  hb2lfyf  +  0(h 3) 

-  Assembly: 

y(t0  +  h)  —  y(tQ )  +  hw \kx  +  hw2k2  +  0(h 3 ) 

y(tQ+h)  =  y(t0)  +  h(\vx+w2)f  +  h2w2a2ft  +h2w2b2lfyf  +  0(h 3) 


-  Constraints  defining  second  order  methods: 

wx  +  w2  =  1  w2a2  -  w2b2l  -  j 
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Second  order  methods  Third  order  method 


0 

O 

a 

a 

1 

2 

1 

2 

i- 1  1 

2  a  2  a 

1 

—  12 

1  2  1 

6  3  6 

An  embedded  second  order 
method  with  a=1/2 


Subtraction  of  two  methods  approx,  error  within  step. 

Error  is  a  state  vector  difference.  „„„ 
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Application 

-  Fixed  step 

-  Adaptive  step 
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Case  1  Simulation  Step  Size 
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Description: 

•  Use  a  fixed  step  integration  method. 

•  Monitor  error  at  each  DoF. 

•  Error  exceeding  tolerance  indicates 
the  system  is  too  stiff  to  integrate 

time  (seconds) 

accurately. 

•  Substitute  impulsive  event  response  in  zero  time. 

•  Continue  integrating. 


•  Determine  accurate  loads  independently  from  initial 
conditions  for  impulsive  event. 
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MODELING  rho  SIMULATION,  TESTING  rkd  VALIDATION 


Generalized  Momentum  and  Impact  Solutions: 


I 


+  dv.  +  dco, 

ni:V: - -  +  1 :  *  CO:  * 


\ 


l  l 


du 


l  l 


du 


r  J 


dv 

-  FcAt - 

du.. 


Z 


-  dv.  T  _  dcO; 

m  :V: - -  +  I:  'CO:  * 


A 


l  l 


i  V 


du 


l  l 


du 


r  J 


+  ^ 

ve  n 


sv'c  ■  h 


14  August  2010  UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


16 


GVSETS 


Ks.t*v 

k  &\ 

r 


Hybrid  Step  Method 


* 


'in 

V 

MODELING  rho SIMULATION,  TESTING  m VALIDATION  '  H 


DESIRED  Loads  Method 


Hybrid  Impulsive  Load  Time  History  External  Integration 


Applied  Finite  Time  Impulses 

•  A  compromise,  in  some  sense  a  more  “general”  solution. 

-  Modeling  the  design  intent  of  a  component  throughout  a  stiff 
region  with  prescribed  motion  constraint. 

-  Peak  loads  reconstructed  using  initial  conditions  for  impulse. 


Hybrid  Step  Method 


Road  Arm  Rates 
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Hybrid  Step  Method 


Road  Arm  Rates 


COR  =  1 
Better  values 
needed 
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•  A  new  method  for  rapid  loads  evaluation  has  been 
presented  in  combination  with  real-time  requirements. 

•  Stiffness  is  identified  dynamically  rather  than 
geometrically  (no  contact  detection). 

•  Zero  time  events  proved  difficult  and  are  the  subject  of 
on  going  work. 

•  Finite  time  impulses 

-  Have  broader/easier  application. 

-  Peak  loads  may  be  post  processed  off-line  and  independently. 
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